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Abstract. Photoevaporation driven by the ultraviolet ra- 
diation from massive stars severely limits the lifetime of 
protoplanetary discs around stars formed within stellar 
clusters. I investigate the resulting influence of clustered 
environments on the probability of giant planet formation, 
and show that for clusters as rich, or richer than, Orion, 
the time available for planet formation is likely to be lim- 
ited to the length of any delay between low mass and 
high mass star formation. Under popular models for the 
formation of massive planets, the fraction of stars with gi- 
ant planets in rich clusters is expected to be substantially 
suppressed as compared to less clustered star formation 
environments. 

Key words: accretion, accretion discs - solar system: for- 
mation - stars: formation - planetary systems - globular 
clusters: general - open clusters and associations: general 
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1. Introduction 

Almost all low mass stars are born in binary systems 



(Reipur th fc Zinnecker |1993| Ghez et al. |l997l ; Kohler & 
Leinert 1998), which are themselves often part of larger 
stellar clusters (Clarke, Bonnell & Hillenbrand 2000). In 
some well-studied star forming regions, such as Taurus- 
Auriga, these are merely loose aggregates of young stars 



(Gomez et al. 1993), which may be expected to have little 
influence on the evolution of protoplanetary discs. This 
is not, however, the case in richer clusters - of which 
Orion is the nearby prototype - where the radiation and 
winds from massive stars have a dramatic impact on the 
local environment (Palla & Stabler 1999; Johnstone & 



Bertoldi 2000; Chevalier 200C). Planet formation in the 



discs around low mass stars in such clusters cannot then 
be considered in isolation. 

Observations of discs in Orion suggest that photoevap- 
oration in the radiation field of the massive stars is likely to 
be the dominant process leading to the destruction of discs 



(Johnstone, HoUenbach & Bally 1998). In this paper, I as 



conclusion is that the formation of giant planets in clusters 
richer than Orion is likely to be heavily suppressed, un- 
less either they form from prompt hydrodynamic collapse 
(e.g. Boss 1998), or there is a substantial delay between 
the onset of low mass and high mass star formation. Pre- 
liminary indications from a Hubble Space Telescope search 
for planetary transits in the globular cluster 47 Tucanae 
(Gilhland et al. |2000| ; Brown et al. pOOOD , in which the 
fraction of giant short-period planets appears to be at 
least an order of magnitude below the value in the so- 
lar neighbourhood (Marcy & Butler 2000D , are consistent 
with this conclusion (provided, of course, that we assume 
that planet-bearing stars in the solar neighbourhood were 
not themselves formed within rich clusters). With hind- 
sight, however, numerous plausible explanations for the 
deficit of planets in this dense, low metallicity system are 
likely to be forthcoming (see e.g. Bonnell et al. 2000). 

2. Disc lifetime in stellar clusters 

2.1. Ultraviolet flux 

Theoretical estimates suggest that there is no reason why 
protoplanetary discs around isolated low mass stars should 
not survive for lengthy periods, especially if the rate of an- 
gular momentum transport (and hence accretion) within 
the disc is reduced due to the low ionization fraction at a 
few AU (Matsumoto & Tajim a |1995| ; Gammie |1996| ; Ar- 
mitage, Livio & Pringle 200C| ). Although estimates of the 
ages of pre-main-sequence stars are subject to significant 
uncertainties (Tout, Livio & Bonnell 1999), some low mass 
stars (M ss 0.5Mq) do appear to retain their discs for 



more than 10 Myr (Strom |1995| ; Brandner et al. |2000| ) . The 
disc lifetime is found to be much shorter for more massive 
stars (Strom |1995 ), and this, coupled with the observa- 
tion of substantial mass loss rates from discs in the Orion 



Nebula (McCuUough et al. 1995; Johnstone, HoUenbach & 
Bally 1998; and references therein), suggests that photo- 
evaporation of the discs is an important process that can 
lead to their destruction (for a review see e.g. HoUenbach, 
Yorke & Johnstone pOOOj) 



sume that this is the case, and discuss the implications for 
planet formation in clusters of varying richness. The basic 



Photoevaporation of discs has been extensively stud- 
ied, both when the source of ionizing radiation is the cen- 
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tral star (Shu, Johnstone & HoUenbach 1993; HoUenbach 



et al. 1994), and for the case relevant here where the radia- 
tion field arises externally (Johnstone, HoUenbach & Bally 



199S; Storzer & HoUenbach 1999). In the simplest analy 



sis, ultraviolet radiation heats the disc surface, raising the 
sound speed to c^. Beyond a radius, 



Rg « 0.5- 



.GAL 



(1) 



the heated gas is unbound and flows away as a ther- 
mal wind. Both Lyman continuum EUV photons (A < 
91.2 nm, hi/ > 13.6 eV), and less energetic far-ultraviolet 
(FUV) radiation (6 eV < hv < 13.6 eV) can contribute 
to disc mass loss. However, for « M©, FUV-driven 
flows (for which the sound speed in the heated layer is 
^ 3kms~^, as compared to ~ lOkms"^ for EUV-driven 
flows) have i?g that is substantially larger than the radii 
of principal interest for planet formation. In this paper 
we therefore concentrate on EUV-driven mass loss, while 
noting that the effects of FUV radiation will also be im- 
portant for lower mass stars. 

To estimate the strength of ultraviolet radiation in 
clusters of varying richness, we first need the emission 
from an individual star. For stars with 0.8 Mq < < 
120 Mq, we obtain the luminosity and effective temper- 
ature Tf, from stellar models computed by Schaller et al. 
( |1992D for a metallicity Z — 0.02 (for these purposes, dif- 
ferences in metallicity are of secondary importance). We 
use the model output closest to t = 10^ yr, which should 
be appropriate for young clusters of any age given the very 
short pre-main-sequence phase of the massive stars that 
dominate the UV flux. The model output is then combined 
with theoretical stellar atmosphere models by Kurucz (e.g. 
Buser & Kurucz |l992|) to yield the EUV and FUV output 
as a function of stellar mass. 

Fig. 1^ shows the EUV and FUV luminosity as a func- 
tion of stellar mass. A negligible fraction of the total flux 
is emitted in the EUV by stars with Te < 2.5 x 10"* K, 
which corresponds to a mass k, 12Mq. The fraction 
of flux in the EUV band then rises roughly linearly with 
increasing Tg to ^ 0.6 of the total for Af* ~ lO^M© stars 
with Te « 5 X 10'' K. Significant fluxes of FUV radiation 
are produced by lower mass stars with > 10* K. 

For a standard mass function the most massive stars 
are rare. Fig. || also shows which stars dominate the output 
of EUV and FUV radiation in a cluster. We assume a 
Kroupa, Tout & Gilmore (199C) mass function for < 
1 Mq, and a Salpeter (|1955D form, n{M)dM oc M^^-^s^ 



for higher masses. Both the FUV and (especially) the EUV 
flux is expected to be dominated by the most massive star 
in a cluster. 

Using this mass function, with a low mass cutoff at 
0.1 Mq, and a high mass cutoff at 120 Mq, we generated 
realizations of clusters with number of members ridustor 
ranging from 10 to 10^. Fig. || shows the resultant in- 
tegrated EUV luminosities. The median luminosity rises 
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Fig. 1. Upper panel, the luminosity in EUV (solid lines) 
and FUV (dashed lines) for stars of mass M*. The lower 
panel shows, for a Salpeter mass function, the expected 
relative contribution to the total flux from stars of differ- 
ent masses. The units in the lower panel are arbitrary. 

steeply for smaller clusters, which do not fully sample the 
high mass end of the mass function, and becomes linear 
only for larger clusters with ^duster > 10"^. The depen- 
dence of the luminosity on the single most massive star 
that a cluster happens to have implies a large dispersion in 
the integrated luminosity. The range of luminosities that 
encompasses 90% of the probability distribution spans al- 
most 3 orders of magnitude at ndustcr = 10^, and is still 
a factor of 3 at riciustcr — 10*. This implies that other- 
wise similar clusters in which the maximum stellar mass 
varies can have widely different UV environments, lead- 
ing to substantial dispersion in the predicted lifetimes of 
circumstellar discs. 

2.2. Disc lifetime 

The theory of EUV-driven flows from discs provides a pre- 
diction for how the mass loss rate scales with the flux of 
ionizing photons. In the simplest analysis, the mass loss 
rate for a disc of fixed radius, at distance d from a source 
emitting a flux of ionizing photons <&, scales as (Bertoldi 
& McKee |1990] Johnstone, HoUenbach & Bally |1998D , 



(2) 



More sophisticated analysis are available (Storzer & Hol- 



lenbach 1999; Richling & Yorke 1998), but are hardly war- 



ranted here given the uncertainties. 

To calibrate the predicted disc lifetime, we make use 
of the results obtained in Orion. The most massive star in 
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Fig. 2. The integrated EUV luminosity for clusters with 
v^ciuster stars. The points show median values, the error 
bars show the range that encloses 90% of the distribution 
of luminosities. 



the Orion Trapezium, 9^ Ori C, has a mass estimated at 
around 40 — 50 Mq (Hillenbrand 1997). Both this mass. 



and the estimated ionizing flux in Orion of the order of 
$ ~ 10"^ s-\ fall withing the range expected for clus- 
ters with n ~ 10'^. Within Orion, one of the best studied 
evaporating discs, HST 182-413, at a projected distance 
of 0.12 pc from 0^ Ori C, has an estimated mass loss rate 
-^outflow = 4.1 X 10"'^ M0yr~-^. The disc radius rdisc is 
~ 100 AU, and the estimated disc mass 0.04 Mq (John- 



stone, Hollenbach & Bally |1998[ ). This disc mass is consis- 
tent with the upper end of the distribution of disc masses 
inferred fr om mm - wavelength observations by Osterloh & 
Beckwith ( 1995 ). For this object, the characteristic disc 
hfetime is therefore idisc = Mnitiai/Moutfiow ~ 10^ yr. 

For planet formation, we are interested in smaller discs 
than that which presently surrounds HST 182-413. These 
will have smaller mass loss rates, which in the case of EUV- 

'^di'sc (Johnstone, Hollenbach & 



Bally 1998). Scaling to a solar system sized disc (30 AU), 
and assuming conservatively that the mass remains similar 
to that of HST 182-413, we obtain tdisc ~ 6 x lO'^ yr as a 
simple estimator of the disc survival time for <i> — 10*^ s~^ 
and d = 0.12 pc. This estimate of the time available for 
planet formation before the disc is destroyed is certainly 
crude, and in particular ignores completely the effects of 
angular momentum transport in the disc. However, viscos- 
ity will tend to hasten disc destruction by moving material 
to larger radii where it can be lost via photoevaporation 
more easily. 



Fig. 3. The predicted dependence of the disc lifetime on 
cluster size, for discs at different radii, 0.1 pc (lower dashed 
line), 0.3 pc (dotted line), 1 pc (long dashed hue), and 3 
pc (dot-dashed line). The error bars plotted for the 0.1 
pc curve show the range of predictions that include 90% 
of the probability distribution. Errors on the other curves 
are identical. 

Fig. ^ shows the dependence of the predicted disc life- 
time on the size of the cluster, for discs located at dis- 
tances d from the principal source of ionizing radiation. 
For smaller clusters this will typically be a single mas- 



sive star, while for rich clusters with ridustcr 



10-* 



larger several stars will contribute. We assume that, as 
in Orion, these stars are formed near the cluster centre 



(Bonnell & Davies 1998), and thus can be treated as a 
single radiation source. If in reality the massive stars were 
instead distributed throughout the cluster, that would re- 
sult in a larger volume over which significant ablation of 
discs would occur. Similarly, mixing of the cluster over 
timescales shorter than the disc lifetime increases the frac- 
tion of stars whose discs are affected by photoevaporation. 

By construction, the dependence shown in Fig. ^ is 
consistent with a destruction timescale of a few xlO^ yr 
in the inner 0.1 pc) part of Orion near 6^ Ori C. It is 
also consistent with observed disc lifetimes of a few Myr, 
extending upwards to more than 10^ yr, in poor clusters 
with just a few hundred members. A large dispersion in 
disc lifetime is expected in this case. 

For clusters richer than Orion, the volume over which 
photoevaporation significantly curtails the disc lifetime 
grows. For large ndustcr, $ oc riciuston and the volume 
within which mass loss is significant scales as 'T-chfstor- 
Thus, for rich clusters, a larger fraction of the discs will be 
ablated by the action of the radiation field, even if the cen- 
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tral cluster density remains constant. For ^duster = 10^, 
the disc lifetime is expected to be severely curtailed (less 
than a few x 10^ yr) out to distances of the order of 1 pc, 
and to be significantly reduced (to < 10^ yr, an order of 
magnitude below the value inferred for isolated discs), out 
to distances of perhaps 3 - 5 pc. 

2. 3. Planet formation 

In the Solar System, the observation that the gas fraction 
of the outer planets varies substantially has been taken 
to imply that these planets formed at roughly the epoch 
when the gas disc was being dissipated, probably at about 
a few xlO^ — 10^ yr (e.g. S hu, J ohnstone & HoUenbach 



star formation appears to have been underway well before 
the formation of the high mass members of the cluster 



(Palla & Stabler 1999), while in some theoretical models 
for high mass star formation prior low mass star forma- 
tion is a necessary ingredient (Bonnell, Bate & Zinnecker 



1998). More generally, unless star formation in a cluster is 



somehow synchronized to better than the sound crossing 
time of the gas in the star forming region, there is bound 
to be some spread in the times at which stars form. Since 
high mass stars are relatively scarce, this would typically 
lead to some low mass stars in the cluster forming well be- 
fore the first massive star turns on and begins the process 
of disc destruction. 



1993| ; Storzer & HoUenbach [19991) . A few Myr is also the 3, Discussion 



typical theoretical estimate for the time required for a 
growing core at 5 AU to reach a mass where it can rapidly 
ac crete g as from the disc in a runaway manner (Pollack et 
al. 1996). Somewhat shorter timescales are derived if the 



cores initially migrate inwards through the gas disc, but 
subsequently halt an d accr ete the gas at smaller radii (Pa- 
paloizou & Terquem 199S| ). These timescales suggest that 
a disc lifetime of 10^ yr might well pose difficulties for gi- 
ant planet formation, at least in some theoretical models, 
while lifetimes closer to 10^ yr would very probably pre- 
clude giant planet formation. From Fig. ^, this suggests 
that giant planet formation will be strongly suppressed in 
clusters with ^ 10^ stars out to around 1 pc, and possi- 
bly at substantially larger radii. We also note that since 
short period planets (analagous to that orbiting 51 Peg) 
are unlikely to have been formed in situ (Bodenheimer, 
Hubickyj & Lissauer 200C), the disc must survive for a 
significant additional time, subsequent to planet forma- 
tion, to allow inward migration through the gas disc (Lin, 
Bodenheimer & Richardson 199(^ ) . Some additional reduc- 
tion in the frequency of such systems would then result. 
Low mass planets are largely unaffected by these consider- 
ations. Indeed, since rapid inward migration through the 
disc can substan tially deplete the population of such ob- 
jects (e.g. Ward 1997 , and references therein), the rapid 
removal of the gas could in principle even enhance the 
initial frequency of low mass planets in clusters. 

The conclusion that rapid disc destruction in clusters 
severely reduces the probability of planet formation can 
be evaded in two ways. First, massive planets might form 
from the disc via direct hydrodynamic collapse (e.g. Boss 



199S, 2000; Cameron 1978). If this process occurs at all, it 



is most likely at very early times when the disc is massive 
and vulnerable to gravitational instability. A detection of 
massive planets around systems where the disc lifetime 
was very short would constitute indirect but persuasive ev- 
idence for the importance of this process. Second, if there 
is a significant delay between the epoch of low mass and 
high mass star formation, planets could have time to form 
before the discs began to be exposed to ionizing radiation. 
In Orion, observations suggest that at least some low mass 



As searches for extrasolar planets extend beyond the im- 
mediate solar neighbourhood, it will become possible to 
study how the frequency and properties of extrasolar plan- 
etary systems depend upon the star formation environ- 
ment. In this paper, I have argued that in what may be 
the typical setting for star formation (Clarke, Bonnell & 
Hillenbrand ^000 ) - a cluster containing a mix of high and 
low mass stars - photoevaporation of protoplanetary discs 
severely limits the time available for giant planet forma- 
tion. Observations indicate that discs are being rapidly 
destroyed by this process within a few tenths of a pc of 
Orion's Trapezium, and the effect will be stronger in richer 
clusters. For a cluster of 10^ stars, the disc lifetime is likely 
to be significantly reduced when compared to that of discs 
around isolated stars, out to a distance from the cluster 
centre of several pc. Unless planets form contemporane- 



ously with the disc itself (Boss 1998), or there is a signif- 
icant delay (perhaps 10^ yr or more) between the epochs 
of low mass and high mass star formation, the fraction 
of stars in rich clusters with massive planets is likely to 
be small. This provides both a possible explanation for 
the apparent dearth of ma ssive, short-period p lanets in 47 
Tucanae (Gilliland et al. |2000| Brown et al. |200C| ), and 
suggests that planet formation could also be noticeably 
suppressed in substantially less dense clustered environ- 
ments. 
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